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Abstract 

I> 

ON 

{S| Plasmodium of acellular slime mould Physarum polycephalum is a very large eukary- 

otic microbe visible to the unaided eye. During its foraging behaviour the Plasmod- 
ium spans sources of nutrients with a network of protoplasmic tubes. In this paper 
we attempt to address the following question: is slime mould capable of computing 
transport networks? By assuming the sources of nutrients are cities and protoplas- 
mic tubes connecting the sources are motorways, how well does the Plasmodium 
approximate existing motorway networks? We take the Netherlands as a case study 
for bio-development of motorways, while it has the most dense motorway network in 
Europe, current demand is rapidly approaching the upper limits of existing capacity. 
We represent twenty major cities with oat flakes, place Plasmodium in Amsterdam 
and record how the Plasmodium spreads between oat flakes via the protoplasmic 
tubes. First we analyse slime-mould-built and man-built transport networks in a 
framework of proximity graphs to investigate if the slime mould is capable of com- 
puting existing networks. We then go on to investigate if the slime mould is able 
calculate or adapt the network through imitating restructuring of the transport 
network as a response to potential localalized flooding of the Netherlands. 
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1 Introduction 



The approximation or computation of shortest path transportation networks 
has drawn significant attention from the field of Unconventional Computing 
Sciences. Nature-inspired computing paradigms and experimental implemen- 
tations have been successfully applied to calculation of a minimal-distance 
path between two given points in a space or a road network. Computa- 
tional models of ant-based optimisation have been shown to be an effective 
way of developing novel approaches towards load-balancing of telecommuni- 
cations [12J, which indeed involves dynamical design of transport links for 
packets. Other work includes a shortest-path problem solved in experimen- 
tal reaction-diffusion chemical systems [lj, gas-discharge analog systems [24], 
spatially extended crystallization systems [7J , formation of fungi mycelian net- 
works [15] and Plasmodium of Physarum polycephalum [20] . 

Amongst all experimental prototypes of path-computing devices slime mould 
Physarum polycephalum is perhaps the most cost efficient biological substrate 
available, coupled with fact that it is both easy to cultivate and observe, it 
makes an excellent computational substrate. These are the main reasons we 
adopt it for this work. 

Acellular slime mould Physarum polycephalum has quite a sophisticated life 
cycle [27], which includes various stages such as: fruit bodied, spores, single-cell 
amoebas, and synt syncytium. Plasmodium is a vegetative stage of Physarum 
polycephalum, it is a syncytium, a single cell, where many nuclei share the 
same cytoplasm. The Plasmodium consumes microscopic particles, and dur- 
ing its foraging behaviour the Plasmodium spans scattered sources of nutrients 
with a network of protoplasmic tubes. The protoplasmic network is usually 
optmized to cover all sources of food while still managing to guarantee robust 
and quick distribution of nutrients in the Plasmodium body. Plasmodium's 
foraging behaviour can be interpreted as computation, with data represented 
by spatial distribution of attract ants and repellents, and results represented by 
the structure of protoplasmic networks [10] . Plasmodium is capable of solving 
computational problems with natural parallelism, namely shortest path [20] 
and hierarchies of planar proximity graphs [5j, computation of plane tessel- 
lations [26], execution of basic logical computing schemes [29, 9j, and natural 
implementation of spatial logic and process algebra [25] . For further examples 
see the overview of Physarum-based computers in [10J. 

In previous work [4] we have evaluated the road-modeling potential of P. 
polycephalum, however, previous results were inconclusive. A step forward 
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biological-approximation, or evaluation, of man-made road networks was done 
in our previous work on approximation of United Kingdom motorways and 
Mexican Federal highways by Plasmodium of Physarum polycephalum [8lfTT] . 
In both cases it was shown that transportation links constructed by Plas- 
modium sufficiently determines man-made motorways, with some pernicious 
differences. Comparing results for United Kingdom and Mexico we found that 
shape of a country and spatial configuration of urban areas or cities sufficiently 
determines behaviour of the Plasmodium. More experiments are necessary to 
provide generalisation, in order to develop a theory of slime-mould based road 
planning and urban development. 

In this paper we hope to move towards a more general understanding of slime 
moulds capability to compute road networks by investigating the roads in 
the Netherlands. The Netherlands presents an excellent case study as it has 
the highest density motorway network in Europe. Moreover, the demand on 
the system is at levels which are reaching current limits, with a total length 
of 132, 397km and usage of 140 x 10 9 people per km per yeaQ Such high- 
occupancy may pose a need for urgent expansion of the transport networks 
and a better understanding of the limitations to that growth. The Netherlands 
is also at risk of significant floodingp~[ 

The remainder of this paper is structured as follows. We delinate the ex- 
perimental method and setup in section [2j Section [3] presents the principal 
experimental results, which are then analysed in a framework of proximity 
graphs in section [4| Restructuring of Physarum- approximated transport links 
for the case of partial flooding of the Netherlands is described in section [5| 
The paper then concludes with a summary of the work and ideas for further 
studies in section |6l 



2 Methods 

All experiements are conducted with Plasmodium of P. polycephalum that is 
cultivated in a plastic container. The Plasmodium are first placed on paper 
kitchen towels, sprinkled with still water and fed with oat flakej" 3 ] The exper- 
iments are conducted in 120 x 120 mm polyestyrene square Petri dishes with 
rounded corners. The Plasmodium will eventually grow on Agar plates, which 
are cut into the shape of the Netherlands. The Agar plates are formed using 
2% agar gel (Select agar, Sigma Aldrich). 



www . aut osnel wegen . nl 

2 http:/ /urbanflood.eu/default.aspx 

3 Asda's Smart Price Porridge Oats 
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(a) (b) 




(c) (d) 

Fig. 1. Experimental basics, (a) Outline map of the Netherlands with twenty one 
sources of nutrients indicated, (b)-(d) Snapshots typical setups: urban areas are 
represented by oat flakes, Plasmodium is inoculated in Amsterdam, the Plasmodium 
spans oat flakes by protoplasmic transport network. 

All experiments consider the twenty one most populous urban areas in the 
Netherlands (Fig. [jji): 

(1) Leeuwarden (3) Den Helder 

(2) Groningen (4) Lelystad 
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(5) 


Zwolle 


(14) 


Dordrecht 


(6) 


Haarlem 


(15) 


Nijmegen 


(7) 


Amsterdam 


(16) 


Hertogenbosch 


(8) 


Utrecht 


(17) 


Breda 


(9) 


Amersfoort 


(18) 


Tilburg 


(10) 


Apeldoon 


(19) 


Middelburg 


(11) 


Enschede 


(20) 


Eindhoven 


(12) 


Den Haag 


(21) 


Maastricht. 


(13) 


Rotterdam 







Further we refer to the urban regions as U. The regions in U are projected onto 
the gel in the following manner: oat flakes are placed in the positions of each 
region (Fig. [J). At the be ginning of each experiment a piece of Plasmodium, 
usually already attached to an oat flake, is placed in Amsterdam (region 7 in 
Fig- [J). 

The Petri dishes with Plasmodium are kept in darkness, at a temperature 
of between 22 and 25 C°, except for short periods of observation and image 
recording. Periodically the dishes are scanned using an Epson Perfection 4490 
scanner. Scanned images of dishes are enhanced to increase readability of 
the image, this is done by increasing saturation and contrast (saturation is 
increased to 55 and contrast to 40). A total of 62 experiments were conducted. 

To ease understading of experimental images we provide complementary bi- 
nary version of each image, where appropriate. In these images each pixel of 
the color image is assigned a black color if red R and green G components of 
its RGB color exceed some specified thresholds, R> Or, G > 9q and the blue 
component B does not exceed some threshold value B < 9 B ; otherwise, the 
pixel is assigned a white color (exact values of the thresholds are indicated in 
the figure captions as = #b)). 



3 Transport links via foraging 

In the following we present experimental results which show the the Plasmod- 
ium is capable of computing, or calculating, the transport links between each 
of the twenty one most populus areas of the Netherlands. 

In a laboratory experiment, illustrated in Fig. [2j the following chain of events 
unfolds (dynamics of colonisation is schematically represented in Fig. [3]). An 
oat flake colonised by Plasmodium was placed on top of the oat flake represent- 
ing Amsterdam. In 12 hours the Plasmodium follows gradients of chemoattrac- 
tants, links Amsterdam with Haarlem, and propagates towards Utrecht and 
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(a) t =12 h 



(b) t =34 h 




(e) t =34 h (f) t =57 h 



Fig. 2. Illustrative example of Plasmodium development on configuration of cities 
represented by oat flakes: (a)-(c) scanned image of experimental Petri dish. Time 
elapsed from inoculation is shown in the sub-figure captions, (d)-(f) binary images, 
6 = (100, 100,100). 
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(a) 



(b) 



Fig. 3. Diagram of colonisation dynamics derived from experiments Fig. [2] (a) and 
Fig. [i] (b): links developed in 12 hours after inoculation are shown by red solid 
arrows, in 34 hours by blue dotted arrows, in 57 hours by green dashed lines, in 
80 hours by dash-dotted lines. Large mesh-patterned arrows indicate migration of 
Plasmodium outside the country. 

Amersfoort, spreading in all directions except north-west (Fig. [2^id). After 
34 hours the Plasmodium colonizes Leeuwarden and Groningen. It develops 
clearly visible protoplasmic tubes, which represent a transport link Amersfoort 

- Lelystad - Leeuwarden - Groningen (Fig. [^e). In the same time interval the 
Plasmodium colonises Apeldoon and start propagations towards Zwolle and 
Enschede (Fig. ^pe). 

After a total of 57 hours the Plasmodium connects Apeldoon with Enschede 
and Zwolle by protoplasmic tubes and colonised south-west part of the country. 
Namelym, the Plasmodium links Haarlem and The Hague and builds a route 
from The Hague to Middelburg and a link Hague - Rotterdam - Dordrecht 

- Breda - Tilburg - Hertogenbosch (Fig. [^f, Fig. |3^i, green dashed lines). At 
the same time the Plasmodium forms a protoplasmic tube directly connecting 
Amsterdam Den Helder, and The Hague with Hertogenbosch, and develops 
the links Hertogenbosch - Nijmegen and Tilburg - Hertogenbosch - Eindhoven 

- Maastricht (Fig.grf). 

We observe that the dynamics of colonisation is non-uniform (Fig. [3^l). The 
Plasmodium does not spread or diffuse in all directions simultaneously but 
rather colonises north-north-west part of the country first and only then ex- 
plores south-south- west. This may be due to the fact that centres of activity 
(biochemical oscillators) form during propagation, and the contractive waves 
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(a) t =12 h 



(b) t =34 h 




(c) t =57 h (d) t =80 h 



Fig. 4. Plasmodium spreads beyond 'dedicated' experimental domain: 
(a)-(d) scanned image of experimental Petri dish. Time elapsed from inoculation 
is shown in the sub-figure captions, (e)-(h) binary images, O = (100, 100, 100). 

evoked by the oscillators that force the protoplasm to move towards the oscil- 
lators. Therefore, if an oscillator is formed in the north part of Plasmodium, 
then propagation in all other directions would be suppressed. 

The Plasmodium of Physarum polycephalum rarely repeats itself in experi- 
mental trials. The overall or average pattern, as we will discuss further in the 
paper, may be the same but a myriad of variations are possible in the course 
of Plasmodium's spatial development. 

Outperforming (spreading out of the dedicated area) and under-performing 
(not colonising the whole area) are typical examples of the varieties in Plas- 
modium behaviour. These two examples are illustrated and discussed below. 

In a substantial number of laboratory experiments, the Plasmodium did not 



8 



(e) t =12 h (f) t =34 h 




( g ) t =57 h (h) t =80 h 

Fig. 4. Continued. 



stop its foraging activity even when all sources of nutrients were occupied and 
the whole agar plate was explored. As shown in Fig. [|a vigorous Plasmodium 
can spread over surrounding Petri dishes, trying to settle on bare plastic. 

In this experiments Plasmodium starts its colonisation in Amsterdam as be- 
fore. It's colinization is more aggresive in this case and it colonises Haarlem, 
Den Haag, Rotterdam, Dordrecht, Breda, Tilburg, Hertogenbosch, Eindhoven 
within the first twelve hours. A pronounced protoplasm transport link is es- 
tablished connecting these cities in a chain (first 12 hours from the moment 
of inoculation, Fig. ^e). 34 hours after inoculation the Plasmodium sprawls 
from Hertogenbosch to Utrecht, Amersfoort and Apeldoon, and then builds a 
transport link Amersfoort-Lelystad-Zwolle-Enschede (Fig. |4]bf ) . 

Protoplasmic tubes connecting Haarlem, Amsterdam, Lelystad with Den Helder 
are grown simultaneously after 57 hours of the experiment. By the same time 
Plasmodium also connects Nijmegen with Apeldoon (Fig. |4j^g). Protoplasmic 
transport links Den Helder - Leeuwarden - Groningen, Rotterdam - Middel- 
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burg and Eindhoven - Maastricht are developed by the 80th hour of Plasmod- 
ium's foraging activity (Fig.[I]ih). A Schematic illustration of the colonisation 
dynamics is shown in Fig. ^p. 

Plasmodium starts to show overperformance after 57 hours of the experiment. 
It sprawls from Den Helder north-westward and from Enschede south-eastward 
onto bare plastic of the experimental container (Fig. ^pg)- The Plasmodium 
does not propagate on the plastic long enough and retracts in few hours (this 
can be seen in Fig. [4]lh). Another sprawling takes place by the 80th hour of 
experimentation, when plasmodium propagates westward of Middelburg and 
sout-eastward of Maasrticht (Fig.[3]lh). See also diagrams of sprawling outside 
the county in Fig. ^jp. Also notice how the plasmodium dynamically changes 
its foraging strategy (Fig. [I]). It first attempts to colonise cities in north-east 
part of the country but then abandons the attempt and move to north-east 
later via the IJsselmeer lake. 

In some experiments the plasmodium never manages to span all cities, and fails 
to collanize some oat flakes. An example is shown in (Fig. [5]), 65 hours after 
inoculation the plasmodium colonises the majority of the Netherlands and 
establishes a network of protoplasmic tubes over most of the cities represented 
by oat flakes (Fig. [SJ^f). Later it goes into a hibernation stage and forms a 
sclerotium. However, at no moment of its development does the plasmodium 
even approach Middelburg. Such situations are rather atypical and did not 
happen often in our experiments. 

As illustrated above, plasmodium is quite an unpredictable creature and the 
patterns formed by its protoplasmic networks in any two experiments rarely 
match each other exactly. Thus we generalise results of our experiments by 
constructing a probabilistic Physarum graph. A Physarum graph is a tuple 
P = (U, E, w), where U is a set of 21 cities, E is a set edges, and w : E —> [0, 1] 
is a probability- weights of edges from E. For every two cities a and b from U 
there is an edge connected a and b if a Plasmodium's protoplasmic link is 
recorded at least in one of k experiments, and the edge (ab) has a probability 
calculated as a ratio of experiments where protoplasmic link (ab) occurred to 
the total number of experiments k. We do not take into account exact configu- 
ration of the protoplasmic tubes but merely their existence, e.g. protoplasmic 
tubes linking Eindhoven with Maastricht always positioned inside the Nether- 
land territory but corresponding edge in Physarum graph represents the tubes 
by straight line. We also consider threshold 9 G [0, 1] Physarum graphs P(0), 
defined as follows: for a, b G U, (ab) G E if w(ab) > 9. 

Threshold Physarum-graphs extracted from 62 laboratory experiments are 
shown in Figs. [5]and[6j The graph becomes planar only for 9 = ^ (Figs. [5^), 
i.e. when edges occurred in over 15% of the experiments. We can therefore 
infer that the Physarum graph is planar. However, with acquiring planarity 
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(a) t =22 h 



(b) t =43 h 




(c) t =65 h (d) t =22 h 




(e) t =43 h (f ) t =65 h 



Fig. 5. Plasmodium does not always span all cities (sources of food): (a)-(d) scanned 
image of experimental Petri dish. Time elapsed from inoculation is shown in the 
sub-figure captions, (e)-(h) binary images, O = (100, 100, 100). 
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(m) e = 1 (n) e = § 

Fig. 6. Configurations of threshold Physarum-gY&ph P(0) for = jjj§ 3 jj^ 5 • • • 3 §| 
Thickness of an edge is proportional to the edge's weight. 
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the graph becomes disconnected: one node, Den Helder city, becomes isolated. 



The Physarum-graphs become acyclic for 6 — || (Fig.|6ji), i.e., when its edges 
appear as protoplasmic tubes in over 41% of the experiments. When the graph 
becomes acyclic it is split into a set of isolated nodes: Den Helder, Leeuwarden, 
Haarlem, Amsterdam, Utrecht, Amersfoort, Enschede, Middelburg, and two 
additional components. One component is a chain of three cities: Lelystad, 
Zwolle and Groningen. The second component is a tree routed in Tilburg. 
The tree has three linear branches: 

• Tilburg - Breda - Dordrecht - Rotterdam - Den Haag 

• Tilburg - Hertogenbosch - Nijmegen - Apeldoon 

• Tilburg - Eindhoven - Maastricht. 

This tree is a characteristic feature of the Physarum-graph and it appears in 
over 60% of experiments. The tree is 'destroyed' when 6 > || (Fig. 6j), then 
only chains remain, which give away isolated nodes with further increase of 
6. Some chains are more stable then others. Thus, the chain Breda - Dor- 
drecht - Rotterdam - Den Haag appears in over 54% of experiments (Fig. |6|). 
While the chain Dordrecht - Rotterdam - Den Haag appears in almost 60% of 
experiments. 

The experiments have now provided a reasonably consistent set of connections 
between the various urban centres in the Netherlands. The next question is 
to asses how well these Physarum graphs approximate the Netherlands mo- 
torway network. A graph H of Dutch motorways is Fig. [7^. We construct the 
motorway graph H as follows. Let U be a set of urban regions/cities, for any 
two regions a and b from U, the nodes a and b are connected by an edge (ab) 
if there is a motorway starting in vicinity of a and passing in vicinity of b and 
not passing in vicinity of any other urban area c G U. If there is a branching 
motorways, which e.g. starts in a goes in the direction of b and c and at some 
point branches towards b and c, we then add two separate edges (ab) and (ac) 
to the graph H. 

The intersection P(9)f)H of Physarum and motorways graphs is shown in 
Fig. [T^cd for # = 0, ^ and j|. A relaxed probabilistic Physarum graph P(0), 
where an edge appears in the graph if it is recorded in at least one experiment, 
matches the motorway graph H almost perfectly. Just three edges of H are 
not presented in P(0) f| H: 

• (Amsterdam, Der Helden), 

• (Zwolle, Apeldoon), 

• (Rotterdam, Dordrecht) (Fig. [7] 



9 = j|is the highest value for which Physarum graph P(0) remains connected 
(Fig. 5]). Graph P(^) RH loses few more edges (presented in P(^) RH: 
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(a)H (b)P(0)DH 




(c)p(DnH (d)p(i)nn 

Fig. 7. Graph H of man-made motorway network is shown in (a). Intersection 
P(#) 9 — 1, jq and j|, of Physarum P and motorways H graphs is shown in 
(bed). 

• (Den Helder, Lelystad), 

• (Den Helder, Haarlem), 

• (Leeuwarden, Lelystad), 

• (Utrecht, Dordrecht), 

• (Utrecht, Nijmegen), 
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• (Breda, Enschede), 

• (Utrecht, Enschede) (Fig.^). 

As soon as 9 reaches value || the graph P(0) f|H becomes separated on iso- 
lated node Den Helder, and there components: 

• chain Enschede - Den Haag -Rotterdam, 

• cycle Nijmegen - Breda - Middelburg with branches Breda - Hertogenbosch 
- Tilburg, Middelburg - Eindhoven and Nijmegen - Dordrecht - Apeldoon - 
Amersfoort, 

• cycle Leeuwarden - Groningen - Zwolle with branches Zwolle - Enschede 
and Zwolle - Lelystad - Amsterdam - Haarlem, Utrecht (Fig. 



4 Comparing motorway and Physarum graphs with proximity graphs 

We conjectured that Plasmodium of Physarum polycephalum constructs planar 
proximity graphs by its protoplasmic network [5j. A protoplasmic network 
constructed in any particular experiment is planar, a generalised Physarum 
graph P may be non-planar. A planar graph consists of nodes which are points 
on a Euclidean plane with edges which are straight segments connecting the 
points. A planar proximity graph is a planar graph where two points are 
connected by an edge if they are close in some sense. Usually a pair of points 
are assigned a certain neighborhood, and the pair are connected by an edge 
if their neighborhood is empty. Relative neighborhood graph [H], Gabriel 
graph [23], /3-skeletons [16] and spanning tree are most known examples of 
proximity graphs. 

For self-consistency we provide a brief definition of the following graphs: 

• RNG: Points a and b are connected by an edge in RNG if no other point 
c is closer to a and b than dist{a ) b) [28J. 

• GG: Points a and b are connected by an edge in GG if a disc with diameter 
dist(a, b) centered in middle of the segment ab is empty [T31I23] . 

• BS(/3): A /3-skeleton, /3 > 1, is a planar proximity undirected graph of an 
Euclidean point set where nodes are connected by an edge if their lune- 
based neighborhood contains no other points of the given set; parameter /3 
determines size and shape of the nodes' neighborhoods [T6] . 

• MST: The Euclidean minimal spanning tree (MST) [22] is a connected 
acyclic graph which has minimal possible sum of edges' lengths. 

The graphs are related as MST C RNG C BS C GG j28U23lfl4] . 

We constructed a relative neighbourhood graph [28j RNG (Fig. [8^), a Gabriel 
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(a) RNG (b) BS(1.5) 




(c) GG (d) MST 



Fig. 8. Proximity graphs constructed on regions U: (a) relative neighbourhood graph 
RNG, (b) /3-skeleton with control parameter 1.5, BS(1.5) (c) Gabriel graph GG, 
(d) Minimum spanning tree MST 

graph [T31I23] GG (Fig. J^), /3-skeleton ((Fig. [8)3) and a minimum spanning 
tree MST (Fig. |8]l) (we rooted MST in Amsterdam) over nodes correspond- 
ing to centres of urban areas. We then calculated intersections of these graphs 
with the Physarum graph P(0) and the motorway graph H, see Fig. [9j 
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(a) PP|RNG (b)Pf|BS(1.5) (c) Pf)GG 




(g)HflGG (h)HflMST 

Fig. 9. Intersection of Physarum graph P(0) (a)-(d) and motorway graph H (e)-(h) 
with proximity relative neighbourhood graph RNG, /3-skeleton B(1.5), Gabriel 
graph GG and minimum spanning tree MST. 

The following is a list of edges of the proximity graphs that are not present in 
Physarum or motorway graphs: 

• PflMST =MST - {(Den Helder, Haarlem), (Zwolle, Apeldoon)} 
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• PR RNG =RNG - {(Den Helder, Haarlem), (Zwolle, Apeldoon)} 

• PnBS(1.5) =BS(1.5)-{(Den Helder, Haarlem), (Zwolle, Apeldoon), (Utrecht, 
Den Haag)} 

• PflGG =GG - {(Den Helder, Haarlem), (Zwolle, Apeldoon), (Utrecht, 
Den Haag), (Den Helder, Amsterdam), (Den Helder, Lelystad)} 

• HflMST =MST - {(Haarlem, Enschede), (Rotterdam, Hertogenbosch)} 

• HflRNG =RNG — {(Haarlem, Enschede), (Rotterdam, Hertogenbosch)} 

• Hf|BS(1.5) =BS(1. 5)— {(Haarlem, Enschede), (Rotterdam, Hertogenbosch)} 

• HflGG =GG — {(Haarlem, Enschede), (Rotterdam, Hertogenbosch), (En- 
schede, Tilburg), (Rotterdam, Nijmegen), (Amersfoort, Dordrecht), (Dor- 
drecht, Tilburg)}. 

The motorway graph closely matches the spanning tree, relative neighbour- 
hood graph and ^-skeleton. Only two edges (Haarlem, Enschede), (Rotter- 
dam, Hertogenbosch), presented in RNG, MST, BS(1.5) do not exist in H. 
The fact that the relative neighbourhood graph is almost a sub-graph of the 
motorway graph indicates intrinsically logical organization of the transport 
networks in the Netherlands. This is because a relative neighbourhood graph 
is commonly considered to be optimal in terms of total edge length and travel 
distance, and is known to be a good approximation of road networks [301I3T] . 

Said that the transport networks in the Netherlands are redundant — from 
slime mould's point of view — because there is a substantial number of edges 
of H not presented in RNG. 

The same can be said about the Physarum graph, because only edges (Den 
Helder, Haarlem) and (Zwolle, Apeldoon) of RNG are not represented by 
protoplasmic tubes. Under-representations of ^-skeleton and GG in between 
P(0) are much more substantial: three and five edges, respectively. 



5 Flooding 



In this section we decribe the experiments which investigate the capability 
of the Plasmodium to adapt the its network during induced flooding of the 
petri-dish. In essence we investigate how the Plasmodium would calculate and 
adapt the transport network, if the Netherlands were to suffer similar flooding 
of its cities and roads. Experiments on flooding were conducted in 12x12 cm 
Petri dishes. A Petri dish was raised by 1-2 cm in its south-east corner and 
partly filled with liquid (distilled water, either pure or coloured with one drop 
of food coloring) (Fig. 10 =l) . In most cases the flooded area included Middel- 
burg on the south-west and Groningen on the north-east. In the central part 
of the country the flooding often reached Zwolle. The exact flooded area var- 
ied between experiemtns due to slight variations in the thickness of agar gel 
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(a) 



(b) 



Fig. 10. Flooding setup: (a) array of Petri dishes during imitated flooding, (b) Flood- 
ing scheme. Flooded part is filled is shaded (green), direction of outward migration 
are shown by arrows, each large arrow represents 20%, while small arrow 10%. 

substrate and minor differences in inclinations of Petri dishes, the flooding 



scheme shown in (Fig. 10) is rather indicative 



Initially Plasmodium reacts to flooding with increased activity. During the 
first few hours of flooding the Plasmodium typically increases its branching 
at the boundary of the flooded area (Fig. [TT^i). Often there are indications of 
indiscriminate increase of foraging, panic foraging. For example, in Fig. [TTfc 
we can see active sprawling of Plasmodium in the areas around Apeldoon, 
Dordrecht and Enschede: no protoplasmic tubes are formed but rather uniform 
sheets of Plasmodium propagate in these areas. Eventually the activity ceases 
and flooded transport links becomes abandoned (Fig. [TTfc). 

In some cases no 'panic' branching occurs, but non-flooded protoplasmic tubes 
become thicker due to increased transport of nutrients and relocation of masses 



of protoplasm from areas affected by flooding (Fig. 12 a). Often only proto- 



plasmic tubes located along flood line are hypertrophied, for example there is 
an increased thickness of tubes along the route Hertogenbosch - Apeldoon - 
Zwolle - Groningen in Fig. [12b. 



With time, water is absorbed by the gel and is sucked under the gel plate 
due to capillary forces, which in turn causes the overall humidity to increase. 
This is associated with a reduced concentration of nutrients, and an increased 
concentration of metabolites ejected in the agar plate, force the Plasmodium to 
abandon the agar plate and migrate beyond. Often the Plasmodium attempts 
to complete evacuation by crawling onto the water surface (for example see 
Fig.pl). 
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(c) t = 46 h 

Fig. 11. Illustration of Plasmodium behaviour during experiment on flooding. Pho- 
tographs are made at different angles. 




Fig. 12. Compensation of transport networks (a) and increase in borderline transport 

(b). 
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(a) (b) 
Fig. 13. Examples of evacuation. 

Outward migration, calculated in nine experiments, are shown in Fig.fTUfc. We 
conclude that if a flooding of the Netherlands were to happen the major impact 
of migration (60%) will be felt by Western Germany, with just a minor impact 
on Northen France (20%) and a very slight impact on Belgium (10%). If the 
Netherlands were flooded the Western Germany accepts the largest impact of 
mass-migration. 



6 Discussion 



In laboratory experiments with plasmodium of Physarum polycephalum we 
discovered that the Physarum protoplasmic network forms a sub-network 
of man-made motorway networks, i.e. every transport link represented by 
Physarum can also be found as a segment of the motorway network. How- 
ever, the converse does not hold. The motorway network is not a sub-graph 
of the Physarum network; there are edges of motorway graph not represented 
by edges Physarum graph. Transport links Amsterdam to Der Helden, Zwolle 
to Apeldoon, and Rotterdam to Dordrecht are never presented in Physarum 
graphs. This may be interpreted as if the man-made motorway network in the 
Netherlands is redundant from Physarum point of view. 

Also, in many experiments Den Helder city remains disconnected from other 
cities. This is either a mishap of the Physarum approach, e.g. plasmodium 
prefers not to enter narrow peninsula of North Holland, or an indication of a 
somewhat inefficient location of Den Helder. 

The most robust component of the Physarum graph, the component which is 
present in the majority of experiments is a tree with three linear branches. 
First branch is Tilburg - Breda - Dordrecht - Rotterdam - Den Haag, sec- 
ond is Tilburg - Hertogenbosch - Nijmegen - Apeldoon, and third is Tilburg 
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- Eindhoven - Maastricht. A possible explanation would be that relative posi- 
tions and distances between the cities in this trees is optimal for Plasmodium 
physiological functioning. Cities are not close enough to be contaminated by 
products of Plasmodium activity but close enough not to put significant strain 
on the pumping of nutrients between distant parts of Plasmodium's body. 

If we prune the Physarum graph by removing edges which occur in less thana 
quarter of experiments and look at the intersection (i.e. the set of edges present 
in both graphs) of this graph with a graph of motorways we find that intersec- 
tion consist of three disconnected components. The first component is a chain 
Enschede - Den Haag - Rotterdam. The second component is a cycle Ni- 
jmegen - Breda - Middelburg with branches Breda - Hertogenbosch - Tilburg, 
Middelburg - Eindhoven and Nijmegen - Dordrecht - Apeldoon - Amersfoort. 
Finally, the third component is a cycle Leeuwarden - Groningen - Zwolle with 
branches Zwolle - Enschede and Zwolle - Lelystad - Amsterdam - Haarlem, 
Utrecht. With respect to the proximity graph the key finding is that relative 
neighbourhood graphs (which is commonly recognised as a best approxima- 
tion of urban streets and transport networks) are almost (apart of two edges) 
subgraphs of the Physarum graph and motorway graph. 

By physically imitating flooding of some parts of the Netherlands we pre- 
dicted that if a real flooding were to occur, the following events will take 
place: substantial increase in traffic on the parts of motorway networks close 
to the boundary between flooded and non-flooded areas, propagation of the 
traffic congestion to all non-flooded parts of the country, complete paraly- 
sis and abandonment of transport network, migration of population from the 
Netherlands to Germany and France, and Belgium. 

The Plasmodium is incredibly similar, esp. in wave-like behaviour, to sub- 
excitable non-linear media [2] and is mainly guided by gradients of chemo- 
attractants [6]. Based on these two facts we could assume the Plasmodium 
must spread from Amsterdam omni-directionally and then start forming branches 
between cities, as e.g. illustrated in Fig. [T4J Such phenomena do not occur, 
rather the Plasmodium behaves more like a 'single-headed' creature, it chooses 
one direction of movement, explores it, then chooses another direction and ex- 
plores it again, see Fig. [3| Such behaviour of the Plasmodium may explain the 
differences between ideal planar proximity graphs and protoplasmic networks 
constructed by the Plasmodium. 

In future work we plan to undertake more experiments on Physarum-based 
imitation of road formation in other European countries, and may be even 
make experiments at a large scale, where Plasmodium grows over the whole 
Europe. We also plan to evaluate the city state of Singapore, this present a 
fairly unique road network that has been carefully planned and constructed 
within the last 50 years. One further objective would be to try and incorporate 
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(a) t = 1 



(b) t = 2 



(c) t = 3 




(g) t = 8 (h) t = 9 

Fig. 14. Snapshots of a growing spanning tree of U rooted in Amsterdam. Vertices 
active at time step t' are shown by black disc at snapshot t = tt 



three dimensional landscape elevation into our laboratory experiments. ' 
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